The Diffusive Galactic GeV/TeV 
Gamma-Ray Background: 
Sources vs. Transport 



Heinrich J. Volk* 



* Max-Planck-Institut fur Kernphysik 
P.O. Box 103980, 69029 Heidelberg, Germany 



Abstract. The diffuse Galactic 7-ray background, as observed with EGRET on 
CGRO, exceeds the model predictions significantly above f GeV. This is particularly 
true for the inner Galaxy. We shall discuss here the contribution of the Galactic Cos- 
mic Ray (GCR) sources, considered as unresolved, and in addition the possibility that 
the transport of the GCRs out of the Galaxy is not uniform over the Galactic disk. In 
both cases the spectrum of the diffuse gamma raysis harder than the GCR spectrum 
in the neighborhood of the Solar system, as observed in situ. The source contribution 
is a necessary and, as it turns out, significant part of the diffuse background, whereas 
the transport effect is one of several conceivable additional causes for the hard diffuse 
7-ray spectrum observed. 



INTRODUCTION 

The observations of the diffuse Galactic 7-ray emission can be described rather 
well by a suitable model for the diffuse interstellar gas, GCR, and photon distribu- 
tions (e.g. Hunter et al. 1997a). However, above 1 GeV the observed average diffuse 
7-ray intensity, foremost in the inner Galaxy, 300° < / < 60°, |6| < 10°, exceeds 
the model prediction significantly. As far as the energetic particles are concerned, 
there are at least two possible explanations for this discrepancy (e.g. Weekes et al. 
1997; Hunter et al. 1997b, and references therein). The high-energy 7-ray excess 
may indicate that the GCR spectrum observed in the local neighborhood is not 
representative of the diffuse CR population in the Galactic disk. An unresolved 
distribution of CR sources is the other possibility. Since the 7-ray emission is the 
product of the energetic particle intensity on the one hand, and of the gas density 
or the photon density, on the other, it is of course possible that deviations from 
the above model assumptions for these latter densities across the Galaxy can also 
lead to changes in the observed energy spectrum of the diffuse gamma rays. We 
shall not discuss such deviations here. We shall rather evaluate the contribution 
of the sources, assumed to be the ensemble of Supernova Remnant (SNR) shells, 



following a recent calculation by Berezhko & Volk (1999). We shall also consider 
the transport of the particles from the same sources out of the Galaxy to naturally 
increase with decreasing Galactic radius (Breitschwerdt et al., 1991). We shall leave 
aside the possibility of new sources of particles, not known in the neighborhood of 
the Solar system. 

GAMMA RAYS FROM THE ENSEMBLE OF SNRS 

Since at best a handful of shell SNRs could be argued to have been detected up to 
now in gamma rays, we shall ignore their discrete contributions and consider the CR 
sources to be spatially averaged over the volume V g = 2.5 x 10 66 cm 3 of the Galactic 
gas disk, with a radius of 10 kpc and a thickness of 240 pc. The corresponding gas 
mass is M g = 4 x 1O 9 M (Dickey & Lockman, 1990). The source input rate in the 
form of energetic particle energy equals v$n 8 E$n, where we take vsn = 1/30 yr, 
Esn = 10 51 erg. The efficiency per SNR is 5 < 1. The total number of localized 
SNRs which still contain their shock accelerated CRs, called here the source CRs 
(SCRs), is given by N SN = ^snTsn, where T SN is their assumed life time, i.e. the 
time until which they can confine the accelerated particles in their interior. Thus 
Nsn is dominated by the population of old SNRs. We estimate T$n — 10 5 yr. 
After the time T SN the SCRs rather quickly become part of the ordinary GCRs 
that presumably occupy a large Galactic residence volume uniformly. 



Acceleration model 

We assume the overall SCR number inside a single SNR to be given by a power 
law spectrum NscR.dE oc e~ 7SCfl dE in energy e in the relativistic range. 

Averaged over the disk volume, the spatial density nsc R (c) of SCRs is given 
n S cR^) = N SC R(e)N SN /V g , with energy density e SCR = N SN 5E SN /V g . In terms 
of escR, w e have 
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for Jscr > 2. The same expressions hold for the GCRs, given tQCR an d Igcr ■ 
For the SCR we may quite possibly have ^scr = 2, and then 
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where e max ~ 10 5 mc 2 is the maximum SCR energy. 



The 7r°-decay prodution rate is given by 

<5 7 (e) = Z^a pp cN g n(e), (4) 

(Drury et al. 1994), which leads to the ratio R = Q^ CR /Q^ CR of the 7-ray pro- 
duction rates due to SCRs and GCRs, given by 

D , , Z* cr N sn 5Esn 
it(e 7 ) = 



Z^ cr (igcr - 2) ln(e maa; /mc 2 )^e G ci? 



x C ( ) , (5) 

VmcV 

where C is the ratio Ng CR /N g JCR ) Ng CR is the mean source gas density, and N g JCR 
denotes the average gas density in the disk. 

With 5 = 0.2, e GC R — 2 x 1CT 12 erg/cm 3 for the relativistic part of the GCRs, 
and 7gor = 2.75 which results in Z^ CR /Z^ CR = 10 (Drury et al. 1994), we obtain 
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for Iscr = 2. 

The total 7-ray spectrum measured from an arbitrary Galactic disk volume is 
then expected to be 

dN 1 dNl CR . A n . 
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where the additional factor 0.4 is introduced to approximately take into account 
the contribution of GCR electron component to the diffuse 7-ray emission at GeV 
energies, and where dN ^ CR is taken from the paper by (e.g. Hunter et al. 1997b). 



"Leaky Box" -type model 

We can derive very similar results from a leaky box-type balance equation 

^GCfi(e) N SCR {e) 
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where V c (e) is the energy-dependent residence volume occupied by GCRs that reach 
the gas disk during their constant mean residence time r c ~ 3 x 10 7 yrs in V c (e). In 
the case of an extended Galactic Halo, V c (e ^> lGeV) ^> V g (Ptuskin et al. 1997). 
Using eq. (4) we can write 

n SCR _ V c Tsn _ T S n ^ 
nccR V g T c Tg 



The GCR residence time in the disk volume 



Tg = r c V g /V c 
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can be derived from the measured grammage x = 14t>/c(e/4.4 GeV) -0 ' 60 g/cm 2 , 
for e > 4.4 GeV, and x = lAv/c g/cm 2 , for e < 4.4 GeV (Engelman et al. 1990). 

At relativistic energies e > mc 2 , the GCR spectrum and the overall SCR spec- 
trum N SC r oc e'^scR are connected by the relation 



Taking j S cr = 2.15, which leads to Z^ CR /Z^ CR = 7.5 (Drury et al. 1994), we 
obtain for e~ > 4.4 GeV: 



(Berezhko & Volk, 1999). 

The question is, of course, whether the SN confinement time Tsn is time depen- 
dent. Probably this dependence is T S N( e ) = ^o( e / 'e m ai)~ 5 where t is the sweep-up 
time when the SNR enters the Sedov phase and the shock speed begins to decrease 
with time. For average ISM parameters t ~ 10 3 yr. 



In Fig. 1 we show the measurements by Hunter et al. (1997a) and our two 
estimates for the total 7-ray emission, from GCRs plus SCRs. They demonstrate 
that the SCR contribution for the acceleration model exceeds the leaky box values 
for all energies. The reason is that for our empirical model the acceleration efficiency 
for the relativistic part of the spectrum is only 5 ~ 0.08. This is probably due to the 
fact that the mean injection efficiency at the SNR shock is lower than the values 
typically assumed for a parallel shock by a factor of a few. We take the lower 
value for the 7-ray emission in Fig.l as the most reliable estimate for the expected 
diffuse 7-ray emission, including the SCRs. Nevertheless the SCR distribution, 
which is about 10 percent at GeV energies, becomes dominant beyond 100 GeV, 
and exceeds the GCR emission at 1 TeV by almost a factor of 10. It would be very 
interesting to detect the diffuse Galactic 7-ray emission at 1 TeV in order to test 
this prediction. 

Until now we have only discussed the 7-ray emission from hadronic SCRs. In 
fact, there are many reasons to assume that electrons are equally well accelerated 
in SNRs, even if their injection into the shock acceleration process is much less well 
understood. The inverse Compton emission by SCR electrons can be comparable 
with the hadronic emission, even though it does not contribute at TeV energies. 
For a more detailed discussion we refer to the paper of Berezhko & Volk (1999). 



ISCR = Igcr ~ 0.6 = 2.15. 
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RESULTS INCLUDING THE SCRS 
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FIGURE 1. The differential diffuse 7-ray energy ffux vs 7-ray energy above 4.4 GeV (cf. 
Berezhko & Volk, 1999). The heavy symbols arc the EGRET measurements, and the dash-dot line 
is the model prediction of Hunter et al. (1997a). The full curve corresponds to our acceleration 
model with "fscR — 2, whereas the dashed curve corresponds to the Leaky Box model. Both 
theoretical curves incorporate energy-dependent loss from the acceleration region. 



TRANSPORT EFFECTS 

The models used to fit the 7-ray data from, say, EGRET assume a GCR energy 
spectrum that is uniform throughout the Galaxy. This tacitly assumes that the 
GCR transport properties leading to the escape from the Galaxy are everywhere 
the same. However that needs not be the case, and in fact is almost certainly not 
true. The dynamical processes leading to GCR escape depend on the strength of 
the regular magnetic field and on its fluctuation characteristics, as well as on the 
CR pressure, and the gravitational field. An example is the formation of Parker 
bubbles which remove the enclosed CRs through their boyant rise into the Halo 
and ultimately into the Intergalactic Medium. Another example which we wish 
to discuss here in some more detail, involves the Galactic Wind which is partly 
driven by the GCRs themselves (e.g. Breitschwerdt et al., 1991, 1993; Zirakashvili 
et al., 1996). In fact, the wind velocity perpendicular to the disk - in z-direction - 
is much larger in the central regions of the Galaxy than at larger radii, through the 
radial variation of the Galactic gravitational field alone (see Fig. 2). This implies 
that for a given particle energy the boundary seperating the dominantly diffusive 
transport perpendicular to the Galactic disk near the disk from the dominantly 
convective transport at greater halo heights moves down in direction to the Galac- 
tic midplane in the inner Galaxy. Since the GCR diffusion coefficient increases 
with energy, the position of this boundary will depend on energy. As shown by 




FIGURE 2. The terminal Galactic Wind velocity, and the base mass flux density, as functions 
of radius in the Galactic disk, cf. Breitschwerdt et al. (1991). All ISM parameters at the base of 
the wind were considered uniform. The radial variation of the Galactic gravitational field alone 
is sufficient to produce this radial gradient 

Ptuskin et al. (1997), the energy spectrum of the GCRs is typically oc E~ 19 in 
the convection region compared to the standard spectrum oc E~ 21 in the diffusive 
confinement region of volume V c discussed in subsecion 2.1. A line of site that in- 
tersects this boundary will therefore receive gamma rays from two regions of very 
different GCR energy spectra, emitting correspondingly harder spectra than does 
the diffusive confinement region alone. Qualitatively this implies a hardening of 
the truly diffuse 7-ray spectrum with Galactic longitude towards the inner Galaxy, 
for given latitude. However, the effect will disappear for high enough energies when 
the convective zone does no more extend into regions of significant gas density. 

Thus, in contrast to the contribution of the sources, this transport effect looses 
importance at high energies. 

It remains to work out this effect quantitatively. But its very existence illustrates 
the interest we should attach to the measurements of the diffuse Galactic 7-ray 
emission over an as wide as possible range of energies. 



CONCLUSIONS 



The foregoing discussion shows that there are at least two mechanisms of basic 
physical interest that contribute to a deviation of the diffuse Galactic 7-ray emission 



spectrum from what would be expected from CR observations in the Solar vicinity. 
The contribution from the SCRs is an inevitable one and is essentially sufficient 
to explain the data at least for the inner Galaxy; it should be part of the 7- 
ray emission model to begin with. Clearly this does not rule out effects from 
potentially existing new populations of CRs, especially electrons, or the influence 
of an increased strength, for instance, of the Interstellar radiation field. This is 
particularly true for high Galactic latitudes. 
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